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Abstract. We present a new plot for representing Raa{4>) data that em- 
phasizes the strong correlation between high-pj^ suppression and its elliptic 
anisotropy. We demonstrate that existing models cannot reproduce the cen- 
trality dependence of this correlation. Modification of a geometric energy loss 
model to include thermal absorption and stimulated emission can match the 
trend of the data, but requires dNg/dy values inconsistent with the observed 
multiplicity. By including a small, outward-normal directed surface impulse 
opposing energy loss, Ap±n, one can account for the centrality dependence of 
the observed Au + Au elliptic quench pattern. We also present predictions for 
Cu + Cu reactions. 
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1. Introduction 

A theoretical model for RHIC mid- to high-p^ RAAi^P) should reproduce both the 
normalization as well as the azimuthal anisotropy of experimental results; its trend 
must follow the data on a V2 vs. Raa diagram. This is actually quite difficult due 
to the anticorrelated nature of Raa and V2 ; previous models either oversuppressed 
Raa or underpredicted [E El El- In Fig- ffl (a) and (b), we combine STAR 
charged hadron Raa{p±) and V2{p±), PHENIX charged hadron Raa{p±) and V2 
centrality, and PHENIX 7r° {p± > 4 GeV) Raa{4') centrality data [(H. We naively 
averaged the STAR and PHENIX Raa{p±) results to approximately match the p± 
bins of their corresponding V2 measurements. We report the Raa and V2 modes of 
the PHENIX 7r° data. The error bars provided are schematic only. 

Hydrodynamics cannot be applied to mid- to high-p_L particles due to the lack 
of equilibrium. Moreover, a naive application would highly oversuppress Raa due 
to the Boltzmann factors. Parton transport theory attempts to extend hydrody- 
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namics' range of applicability to higher transverse momenta. The Molnar parton 
cascade (MPC) succeeded in describing the low- and intermediate-p_L V2 results of 
RHIC by taking the parton elastic cross sections to be extreme, at ~ 45 mb One 
sees in Fig.n(a) that for the MPC, in this instance run at approximately 30% cen- 
trality, no single value of the controlling free parameter, the opacity, x — J dzatPg, 
simultaneously matches the experimental Raa and V2- 

pQCD becomes valid for moderate and higher p± partons, and models based 
on pQCD calculations of radiative energy loss have had success in reproducing the 
experimental RAAip±) data [ 6 . These models use a single, representative path- 
length; as such, they give V2 ^ 0. To investigate the V2 generated by including 
pathlength fluctuations, we use a purely geometric (neglecting gluon number fluc- 
tuations) radiative energy loss model (GREL) based on the first order in opacity 
(FOO) radiative energy loss equation [U^; it has been shown that including the 
second and third order in opacity terms has little effect on the total energy loss [ 

The asymptotic approximation of this equation is /S.E^^^/ E cc {dNg/dy)L'^ 
We thus use an energy loss scheme similar to [Ej: e — AErad/E = nl. k is a free 
parameter encapsulating the E dependence, etc. of the FOO expression and the pro- 
portionality constant between dNg/dy and ppart- I represents the integral through 
the ID Bjorken expanding medium, taken to be / = /o°°d^ ^j^p-paxtixo+fil)^ where 
Iq = .2 fm is the formation time. We consider only ID expansion here because [ 
IIU) showed that including the transverse expansion of the medium has a negligible 
effect. 

The power law spectrum for partonic production allows the use of the mo- 
mentum Jacobian (p;^ = (1 — e)p^) as the survival probability of hard partons. We 
distribute partons in the overlap region according to pcoW = Taa and isotropically in 
azimuth; hence i?^^((/.; b) = Jd^d;, T^A.,y^b)ii-e{..,y,4';b)r ^ ^^^^^ 4 < n < 5. The 

difference from using n = 4 as opposed to n = 5 is less than 10%, and in this paper 
we will always use the former value. We evaluate Raa{4') at 24 values of <j) from 0-27r 
and then find the Fourier modes Raa and V2 of this distribution. Another method 
for finding V2, not used here, assumes the final parton distribution is given exactly 
by Raa and V2, and then determines V2 from the ratio i?^^(0)/i?^yi(7r/2); this 
systematically enhances V2, especially at large centralities. A hard sphere geometry 
is used for all our models, with Rhs = 6.78 fm ensuring < r"^ >=< >. 
Fig. n (a) shows that even with the HS-geometry-enhanced V2 , the GREL cannot 
recreate both Raa and V2 with a single parameter value. 

2. Exclusion of Detailed Balance and Success of the Punch 

In [12, Wang and Wang derived the first order in opacity formula for stimu- 
lated emission and thermal absorption associated with the multiple scattering of 
a propagating parton, and found AE^f^^/E oc {dNg/dy)L. To model this we use 
g _ AE^^a _ AE^ts — ^/ _ /j/^^ where kI is the same as in the GREL model, 
k is a free parameter encapsulating the proportionality constants in the absorp- 
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tion formula, and I2 represents an integral through the ID expanding medium: 
I2 = J^dl pparti^o + nl). I2 has one less power of I in the integrand; this 
permits a unique determination of the two free parameters, k — .5 and k — .25 
fits the 20-30% centrality PHENIX 7r° RAA{(t>) data point, and allows the model 
to duplicate the data as seen in Fig.n(b). Taking the AE/E equations seriously, 
we invert them and solve for dNg/dy; thus dNg°-'^/dy ~ ^ gTrC^^iii io+^L ^pa^rt , and 

dNg''^/dy ^ ^ z-nCRa'^v2 I^+T-^part; where vi and V2 correspond to the bracketed 
terms in the energy loss and energy gain approximations of [ 9 . For our fit- 
ted values of k, and fc, the choice of £' = 6 GeV, L = 5 fm, and — .4 gives 
dNl'"^/dy ^ 1000 and dNf^/dy ^ 3000 for most central collisions. For £; = 10 
GeV, dN^'"^/dy - 1000 and dNf^dy ~ 9000. The huge increase of dNf^/dy to 
values too large to fit the RHIC entropy data reflects the E^ dependence of the 
Detailed Balance absorption. It seems the only way to have a large enough energy 
gain while maintaining dN^^^ /dy ^ 1000 is to increase as above 1. Note that these 
calculations were performed using a hard sphere nuclear geometry profile, which 
naturally enhances the produced V2 [El- 




Fig. 1. (a) STAR /i± data for 0-5%, 10-20%, 20-30%,. . . , and 40-60%, PHENIX /i± data for 
0-20%, 20-40%, and 40-60%, and PHENIX tt" data for 10-20%, 20-30%,. . . , 50-60% centralities. 
Inability of previous models to fit the data, (b) Addition of thermal absorption or momentum 
punch to GREL; both fit the data, but absorption requires entropy- violatingly large dNg^'^ /dy. 
(c) Cu + Cu predictions for the three models. 

Building on the success of radiative energy loss in reproducing RAAip±), and 
supposing that latent heat, the bag constant, the screening mass, or other decon- 
finement effects might provide a small 1 GeV) momentum boost to partons in 
the direction normal to the surface of emission, we created a new model based on 
the GREL model that includes a momentum "punch," Ap±. After propagating to 
the edge of the medium with GREL, the parton's final, "punched-up" momentum 
and angle of emission are recomputed, giving a new probability of escape. Fitting to 
a single {RaAtV2) point provides a unique specification of k and punch magnitude. 
The results are astounding: one sees from Fig. (b) that a tiny, .5 GeV, punch 
on a 10 GeV parton reproduces the data quite well over all centralities. Fitting 
the PHENIX 20-30% tt° data sets k = .18 and the aforementioned Apj_ = .5 GeV. 
The size of the representative parton's initial momentum is on the high side for 
the displayed RHIC data; however, the important quantity is the ratio Ap±/E. 
Moreover, although the geometry used naturally enhances the V2, we feel confident 
that when this model is implemented for a Woods-Saxon geometry, the necessarily 
larger final punch magnitude will still be relatively small. We expect the magnitude 
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of this deconfinement-caused momentum boost to be independent of the parton's 
momentum; hence V2{p±) wih decrease hke l/p±. Moreover, since e is larger out 
of plane than in, a fixed Ap± enhances RAAi''^/'^) more than RAAiO). These are 
precisely the preliminary trends shown by PHENIX at QM2005. Keeping the same 
values for k, k, Ap±, etc. as for Au + Au, we show in Fig.n(c) the centrality-binned 
Raa and V2 results for Cu + Cu in the three geometric energy loss models. 

3. Conclusions 

By failing to simultaneously match the Raa and V2 values seen at RHIC we dis- 
counted the MPC and pure GREL models. We showed that while including medium- 
induced absorption reproduces the Raa{4>) phenomena, it does so at the expense 
of inconsistent and huge dNg/dy. But the addition of a mere 5% punch created a 
RHIC-following trend. This impulse is small enough to be caused by deconfinement 
effects and future calculations should follow the p±^ dependence of Raa{4>\ P±)- 
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